load eq_SS.mat
idx = find(glob.sf(:,1) == glob.sf(15,1))

subplot(1,3,1)
plot(eq_SS.pf(idx).*(exp(glob.sf(idx,2))))

subplot(1,3,2)
plot(eq_SS.l(idx))

